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Abstract 

We show that the following geometric properties of the motion of discrete and contin- 
uous curves select integrable dynamics: i) the motion of the curve takes place in the 
N dimensional sphere of radius i?, ii) the curve does not stretch during the motion, 
iii) the equations of the dynamics do not depend explicitly on the radius of the sphere. 
Well known examples of integrable evolution equations, like the nonlinear Schrodinger 
and the sine-Gordon equations, as well as their discrete analogues, are derived in this 
general framework. 

1 A historical introduction 

One of the classical problems of the XlX-century geometers was the study of the connection 
between differential geometry of submanifolds and nonlinear (integrable) PDE's. For 
instance, Liouville found the general solution of the equation (known now as the Liouville 
equation) which describes minimal surfaces in E 3 |TJ. Bianchi solved the general Goursat 
problem for the sine-Gordon (SG) equation which encodes the whole geometry of the 
pseudospherical surfaces. Moreover the method of construction of a new pseudospherical 
surface from a given one, proposed by Bianchi ||, gives rise to the Backhand transformation 
for the SG equation Q. 

The connection between geometry and integrable PDE's became even deeper when 
Hasimoto || found the transformation between the equations governing the curvature 
and torsion of a nonstretching thin vortex filament moving in an incompressible fluid 
and the NLS equation. Several authors, including Lamb Q, Lakshmanan Q, Sasaki ||, 
Chern and Tenenblat |]] related the Zakharov-Shabat(ZS) ]l0[ spectral problem and the 
associated Ablowitz-Kaup-Newell-Segur(AKNS) hierarchy |ll]] to the motion of curves in 
E s or to the pseudospherical surfaces and certain foliations on them. 

Almost at that time Sym introduced the soliton surfaces approach, in which the pow- 
erful tools of the 1ST method are used to construct explicit formulas for the immersions 
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of one-parameter families (labeled by the spectral parameter) of surfaces corresponding to 



given solutions of integrable PDE's fl2|| ; see also the recent developments of Bobenko 12]. 



More recently Langer and Perline [14| showed that the dynamics of a nonstretching 
vortex filament in R 3 gives rise, through the Hasimoto transformation, to the recursion 
operator of the NLS hierarchy. Similarly, Goldstein and Petrich [15] showed that the 
dynamics of a nonstretching string on the plane produces the recursion operator of the 
mKdV hierarchy. 

Connections between geometry and integrable PDE's in multidimensions can also be 



found, for example in works of Tenenblat and Terng [16] and Konopelchenko |T7]]. Also at 
a discrete level there is a similar situation. For instance, discrete pseudospherical surfaces 
and the discrete analogues of constant mean curvature surfaces are described by integrable 



discrete analogues of the sine-Gordon (SG) (T^] (T£j] and sinh-Gordon equations [£0[ . Such 
discretizations were found by adapting the main geometric properties of the continuous 
surfaces to a discrete level. 



In two recent papers pij| |22| we have proposed a new geometric characterization of the 
integrable dynamics of a discrete or continuous curve (where, by discrete curve, we mean 
just a sequence of points), based on the following three properties: 

Property 1. The motion of the curve takes place in the iV-dimensional sphere of radius R, 
denoted by S N (R), N > 1. 

Property 2. The curve does not stretch during the motion. 

Property 3. The equations of the dynamics of the curve do not depend explicitly on the 
radius R. 

We remark that Properties 1-3 not only select integrable PDE's, but also provide their 
integrability scheme; in other words, in the process of deriving the dynamics selected 
by Properties 1-3 one discovers "for free" the integrable nature of such dynamics! In 
particular, the spectral problem is given by the Frenet equations of the curve and is a 
consequence of Property 1, and the spectral parameter is given by the inverse of the radius 
of the sphere. 

We also remark that our aproach explains in a simple way Sym's formula [12], which 



allows to calculate, from the wave function of the spectral problem, the surface generated 
by the motion of the curve. 

In papers [^] and (2^] we have dealt with the integrable dynamics of a continuous 
and discrete curve respectively, obtaining, in the case of N = 3, the AKNS [|ll]] and the 
Ablowitz-Ladik (AL) (2*|] hierarchies. As we shall see in the following if, in Property 3, 
we consider discrete time dynamics, one also generates integrable fully discrete evolu- 
tion equations, like the Hirota equation p3|]pq]. Since integrable discrete dynamics can 
always be interpreted as Backlund Transformations (BT's) of the corresponding continu- 



ous dynamics [26], the hierarchies of BT's of integrable systems are also characterized by 
Properties 1-3. 



2 The curve in S N (R) and the associated spectral problem 

2.1 Frenet basis along the discrete curve 

Let us consider a sequence Z 3 ^ r(k) € S N (R) C RA +1 of points of the N-dimensional 
sphere of radius R. We are interested in the sequence r(fc) which gives rise to a piecewise 
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linear curve in S (R). Our goal is to construct an analog of the Frenet basis along the 
discrete curve and of the corresponding Frenet equations. 

By F 1 (k) we denote the 1-dimensional oriented vector subspace of R w+1 given by 
r(k) and let F l+l (k) = F l (k) + F l {k + I). For the point r(k) in general position (what 
we assume in the sequel for simplicity) F l (k) is /-dimensional oriented subspace of R Ar+1 
(I < N + 1). It is also convenient to denote F°(k) = {0}. Now we define the orthonormal 
Frenet basis {^i(k)}iL i n the point r(k) of the discrete curve: fi(k) is the unit vector of 
F l+l (k) orthogonal to F l (k) and correctly oriented. 

The distance A(/c) between points r(k) and r(k + 1) of the curve is given in terms of 
the radius R of the sphere and the angle (fo(k) between fi(fe) and fi(fe + 1) as 

A(k) = RMk) ■ (1) 

We are going to construct N — 1 other angles which play similar role as curvatures of the 
(continuous) curve. 

Both F 2 {k) and F 2 (k + 1) are subspaces of F 3 (k) and their intersection is F 1 (k + l). Its 
orthogonal complement in F 3 (k) is the plane ir\(k). Since ^{k) G F 3 (k) and i2(k)-LF 2 (k), 
then (2(k) € ^i(k); similarly fi(/c + l) € ni(k). By fi(fe) we denote the unit vector of ir\{k) 
normal to ^{k); it is also the vector of vro(fc) = F 2 {k) orthogonal to fo(/c + 1). The angle 
fi(k) between the hyperplanes F 2 (k) and F 2 (k + l) in F 3 (k) (equivalently, between ii(k) 
and f\(k + 1)) is the angle of geodesic curvature. 

Uk+1) 




In general we consider the subspaces F (k) and F (k + 1) of F +1 (k), their intersection 
F l ^ 1 (k+\) and its orthogonal complement 7Tj_i(A;). Since fj(fe) G F l+1 (k) and ii(k)±F l (k), 
then f,(jfe) G 7Tj_i(ife). Similarly, since fj.^jfe + 1) G F l (k + 1) C F i+1 (A;) and fj.^fc + 
ljXF'-^fc + 1), then fj.^jfe + 1) G 7T|_i(A;). 

By fj-i(fc) we denote the unit vector of 7T/_ 1 (/c) orthogonal to f/(/c). One can show 
that f/_i(fc) is also the unit vector of iri^ik) orthogonal to fj_2(fe + 1). This is the 
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consequence of two facts: fj_i(&) G F'(fc) (as fj_i(&) € F l+1 (k) and fi_i(/c)_Lf}(fc)) and 
fj_i(fc)-LF i-1 (A; + 1). The angle </>j_i(&) between the hyperplanes F'(fc) and F l (k + 1) 
of i ? ' +1 (A;) (equivalently, between their normals fi(k) and ?/(&), or between f;_i(/c) and 
f*z-i(^ + 1)) is the angle of the (I — l)th curvature. 

F w (k+1) 




Fig.2 

The transition from the Frenet basis {{i(k)}fL in the point r(k) of the discrete curve to 
{§(k + 1)}j=o * s obtained by the superposition of N rotations of the angles <pi(k) in the 
planes m(k), (I = 0, N — 1): 

{Uk + l)Mk),Hk),-.-,fN(k)} Mk ^ l(k) ... 

... + l),...,f { _ 1 (fc),fK* ; ) J f, +1 (A ; ),... J f JV (A;)} W - l( ^- l( * ) (2) 

{f (fc + 1), H_ a (fc + 1), f,(fc), fiv(fc)} Mk M (k) ... 

... {f (A; + 1), (A; + 1), fv(A:) =f N {k + l)} . 

2.2 The spinor representation of the Prenet equations 

The linear transformation related to the resulting rotation gives the discrete analog of the 
Frenet equations. We will present it using the language of the Clifford algebras and Spin 
groups [27\. Let £ = {ei}fL be a fixed orthonormal basis of R w+1 cosidered as a subspace 



of the Clifford algebra Cl(iV + 1), then any other orthonormal (correctly oriented) basis 
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T = {fi}i =0 can be obtained from £ using an element 5 of the corresponding Spin(A r + 1) 
group 

F = S~ 1 £S . (3) 

Moreover, when another basis T is obtained from T by rotation in the plane (fj, fj) of the 
angle 92, then 

S = (cos | + e iej sin |)5 = 0£S . (4) 

If <S(/c) G Spin(A r + 1) represents the rotation to Frenet basis in point r(k) then it is 
subjected to the equation 

S(k + 1) = Offjiffl • . . . • 0^ {k) 0^ {k) S{k) . (5) 
The arc-length along the curve 

fc— l 

* = E A « ( 6 ) 

i 

in the continuous limit (fo(i) — ► is the arc- length parameter. Moreover 

dS{s) v 5(A; + 1) - / 1 , , \ N 

— lim -i ; = ~B E oi + «i(a)-gi2 + - + k n ^s)E N - 1jN S{s) , 

(7) 

where 

Kl{s)= Hm (8) 

is the /-th curvature of the corresponding continuous curve and Eij = Giej/2 are elements 
of the canonical basis of the orthogonal Lie algebra so(iV + 1) in the Clifford algebra 
representation. The equation (^) is nothing but the classical Frenet equation for the curve 
in S N {R). 

Remark: Even when point r(k) is not in general position, the corresponding spaces F l (k) 
can be defined in a way that their dimension is I. We just keep the space F l (k — p) (p > 0) 
of the nearest point in which it was "properly" defined. In this way one can define, for 
example, the Frenet frame along the geodesic line (which is any big circle): only fo and 
fi vary along the curve, and the rest of the Frenet frame remains as it was defined in a 
starting point. 



3 The discrete curve in S d {R) 

In this Section we investigate in detail the discrete curve in S^{R) with constant distance 
A between the subsequent points. Throught the Section we use the following definitions: 
is = A/ R = ipo(k), ip(k) = (fi(k), 9{k) = 922(^)5 the Frenet basis {^i(k)}f =0 is denoted by 
{r(k),t(k),n(k),h(k)} and consists of the radial, tangent, normal and binormal vectors. 

We show that the Frenet equation (0) reduces to the Ablowitz-Ladik spectral problem, 
and its continuous limit to the Zakharov-Shabat spectral problem. We also present the 
geometric explanation of Sym's formula. 
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3.1 The Hasimoto transformation for the discrete curve in S 3 and the 
Ablowitz-Ladik spectral problem 

It is convenient to modify Frenet basis by a rotation in the normal plane (n(k),h(k)) of 
the angle tr(k) =Ylt~ X (i). 

N(fc) = cos a(k)n(k) - sin a(k)b(k) (9) 

Nj(fc) = sin a(k)n(k) + cos a(k)b(k) . 

This change of basis corresponds to a partial "integration" of the Frenet equations in the 
normal plane, since the vectors N(/c),Nj(fc) does not vary from the point of view of the 
normal plane (this is the discrete analog of the parallel transport in the normal bundle) . 
It is also convenient to interprete any vector of the normal plane as a complex number 

${k) = Re0(&) N (*O + W(fc)Nj(fc) 0(fc) = Re<f>(k) + i Im^(jb) . (10) 

If we define S(k) € Spin(4) by the relation 

H(k) = {r(fc),t(fc),N(fc),Nj(fc)} = S(k)~ 1 £S(k) , 

then 

S(k + 1) = O^Of^O^O^Sik) . 

To represent the above rotation in terms of matrices we first choose the following repre- 
sentation of the basis £ as 4 x 4 Dirac matrices 







(! 


i) 





— «(Tl 


io"i 






where I is the 2x2 identity matrix and 07 are the standard Pauli matrices . The linear 
problem related to the discrete curve takes the form 

a (k + 1) _( S'(k+1) \_(A'(k) \(S'(k) \ 

^ l + j ~l S"(k + 1) ~{ A"{k) J I S"(jfe) i ' 1 j 



with 

K ' I -e i!/ / 2 sin(v9(A;)/2)e- if7 ( fe ) e"^ 2 cos(^(fe)/2) 



e iiy / 2 cos(<p(fc) /2) e"^/ 2 sin(p(fc) /2)e iCT W 



1 / C 



(13) 



wher 



q(k) = t&n(v(k)/2)e ia W , ( = e iu / 2 , (14) 

and A"(k,C) = ^'(^iC -1 )- This linear problem is equivalent |Q| to the Ablowitz-Ladik 
spectral problem pSfl . 

In the limit of the continuous curve we obtain the Zakharov-Shabat spectral problem 

dS'(s) 1 / iX q(s) X 



ds 2\ -q(s) -iX J 
where q(s) = K(s)e ia ^\ a{s) = J s T{s')ds' and A = IT 1 



S'(s) , (15) 
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3.2 The geometric interpretation of Sym's formula 

In short notation: £ «-» {I, ia^ , iu\ , —202} 

H(k) = S"(k)~ 1 £S , (k) 



(16) 



and, as a consequence, for a function q(k), the radius vector of the corresponding curve in 
sphere S 3 (R) is represented by 



r(fc) = RS"(k)- L S'(k) 



(17) 



Suppose one is interested in the radius vector of the curve in R 3 corresponding to q(k). 
One can consider R 3 as sphere of infinit radius but one cannot just take the limit R — > oo 
in the formula above. This way the center of the sphere is fixed while R 3 is pushed away 
to infinity. 




Fig. 3 



To remove this incovenience we first have to shift the basis of R 4 to a point of the sphere: 
r new (k) = r(k) - r . 

If the linear problem is solved under the initial condition S"(0, A) = S"(0,X) = I, 
then S"(k, A) = S'(k, —A). Choosing ro = r(0, A) = Reo = il, one obtains the following 
formula for the Cartesian coordinates (X*(&;)) 3 =1 of the points f (k) of the curve in R 3 



m = j2x*(k) ei = iim l (s'(k, xy's'ikA) -i) = 2s'(k,o)- i ds ' { d k x x) 

i=l " /J ' ' 



A=o 



(18) 



The above formula was first used by Sym [12] in his approach of soliton surfaces. He also 
found its generalization to soliton equations related to linear problems in semi-simple Lie 
algebras. Our approach is more related to the Clifford algebras. 



Remark: A modification of this formula was recently found by Cieslihski [28] in the 
context of conformal geometry of isothermic surfaces in R 3 . 
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4 The integrable dynamics in S S (R) 

In this Section we consider the motion of the discrete curve subjected to Properties 1-3. 
For convenience, we use the following short-hand notation: / for f(k),k € Z and f n for 
f(k + n), n = ±1,±2,... 

The motion of the curve is governed by velocity field v which is convenient to write in 
the form 

sin(AA) 
r, t = v = — 

and, consequently, 



(Vt + Re#N + Im^Nj) , (19) 



S[ t = TS'= l -i^ 1 ^jS' , Tgsu(2) . (20) 

The compatibility condition between equations (12)(|l9|) and ( p0| ) specifies the entries 7 
and 5 in terms of the velocity field: 

6 = iC^-i^+qiV + V{)) 

7 = W + sin(AA)T/ (21) 
Wx-W = -Re(iq(<t> 2 -<f> + qi (yi + V2))) 

and yields the kinematics 

2q yt = -2 cos(AA) (fa + q V l )+TZ (fa + qV{) , (22) 

(l-\q 2 \)V 1 -(l + \q\ 2 )V = qfa+qfa , (23) 

where 

Kf := (l+\q\ 2 ) (/1 + /-1 + 2( Ql E - q^)(E - l^Re-^L^j -2ig(E-l)- 1 Im(g 1 /-g/ 1 ) 

^ " (24) 

and £7 is the shift operator along the discrete curve: Ef = f\. 
Substituting the ansatz 

( \ J=£(cos(AA))^ ^ J (25) 

into equation ( p2| ) and requiring independence of cos(AA), we finally obtain the following 
class of integrable dynamics: 

q, t = ho(n)(l + \q\ 2 )(q 1 -q. 1 ) + h 1 (TZ)(iq) , (26) 

where ho and hi are arbitrary entire functions with real coefficients. 

We remark that equation (^) implies the following interesting connection: 

=fa^ ) +qVi m) (27) 
between the integrable commuting flows 

=K m - 1 (l + \q\ 2 )(qi-q^i) and/or K^=TZ m (iq) , m>0 (28) 
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and the velocity fields. In the continuous limit this reduces to the result of Langer and 
Perline @. 

The simplest examples are the following: 

i) If ho = 1, hi = 0, then 

sin(AA) . _ . sin(AA) . . , 

v= \ J (t-q^)= K x ; (t-|g-i|n) , (29) 

q,t = (I + \q\ 2 )(qi ~ q-i) =: . (30) 

ii) If ho(x) = x, h\ = 0, then 

v = 4sm ^ AA) (Ycos(AA) + ^ (gg_i + qq-x)\ t + A , (31) 

-g_i cos(AA) + ~ ((1 + k-i| 2 )(<7 - q-2) ~ + 99-l))^ 

9,t = (1 + |9| 2 ) ((1 + l9i| 2 )92 - (1 + |9i| 2 )9-2 + 9(9? " 9 2 1) + g(«lg-l - 9-l9l)) =: ^ {2) • 

(32) 

iii) If ho = and h\{x) = x, we obtain 

2sin(AA) 2sin(AA) ^ 
v = (zg_i) = tan( — )b , (33) 

q,t = i(l + \q\ 2 )(qi + q-l) ■ (34) 

This equation has also recently appeared in conection with the Heisenberg XXO antifer- 
romagnet model [29]. 

In the continuous limit, 1Z — 2 reduces to the recursion operator of the continuous NLS 
hierarchy. Moreover the following combination of equation (|34|) with the "zero order flow" 
q,t = iq : 

q,t = i[qx-2q + q-\ + |g[ 2 (9i + 9-1)) (35) 



reduces [23] in the continuous limit, to the NLS equation 

iq,f = q, ss + \\q\ 2 q > t' = -A 2 t . (36) 

which describes the motion of a vortex filament in the localized induction approximation 
|||3^]. We remark that, in this approximation, the velocity field which governs the motion 
of the vortex depends on its curvature k through the relation 

r t = Kh ; (37) 

therefore, since equation (|35|) has in R 3 a velocity field of the same type: 

r t = 2Atan(^l)b , (38) 

we expect it to be a good candidate for describing the motion of a discrete vortex in the 
same approximation. 

Consequently, the continuous limit of the linear combination 

q t = - 2K^ (39) 
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reduces to the complex mKdV equation 

q,t> = q,sss + ^\q\ 2 q,s , t' = 2A 3 t . (40) 

We remark that, in the degenerate case of the curve on 5 2 , 9 = and consequently 
q, (j) £ R. In this case we are forced to choose h\ = and only the first hierarchy survives. 

Property 3 can also be satisfied through a mechanism (different from that of equation 
(|25|)) which gives integrable dynamics with sources (see [ ^H [p2^ j ) . 

5 Discrete-time dynamics and the Hirota equation 

In this section we consider the discrete curve on S 2 (R) with constant distance A = vR 
between subsequent points "moving" in discrete time. 

Now we use the following notation: / for f(k, I), G Z, and iorf(k + m,l + n). 
The Prenet equations read 

5i = Of 2 O^S = AS . (41) 

The discrete-time kinematics can be described in terms of angles to and [i which are 
analogues of the direction and the length of the velocity field respectively. 




The induced kinematics of the Frenet frame reads 

S 1 = O^O^S = BS , (42) 

where the angle p should be calculated from the compatibility condition 

A 1 B = B\A . (43) 

It turns out that it is convenient to use, instead of p, another angle tp = uj 1 + p, which is 
the angle of curvature of the curve {r n } ng z. 
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The compatibility condition ( (43[) written in terms of angles splits into three equations 

ui ui + ip— (uj{ + If 1 - H uj - (uj 1 -ip) 
cos — cos = cos — cos , 

. Ux LOi + (f + U}\ + (fi 1 - lb\ U 0J + UJ l -l{) 

sm — cos = cos — cos , 44) 

2 2 2 2' y ' 

. Hi . UJi + If + Cjj + ip 1 - Ipl , . jUl . CJi + V? - (cjj + - ipi) 

sm — sm i h « cos — sm ± = 

2 2 2 2 

- iu ( ■ V ■ u + ui l -ip . h • ^-(u^-VA 

e sm — sm h i cos — sm 

^ 2 2 2 2 j 

In this paper we consider only the motion subjected to the condition \i =const. This is 
(together with the previous condition v =const.) the discrete analog of the Tchebyschev 
net condition, which in the continuous case gives the sine-Gordon equation |3lf| . The first 
two equations of (PHI) imply 



^ 1 -i-</C 1 = -(^r 1 -^Ci) (45) 
which asserts the existence of the "potential" </>: 

if = <p 2 - 4> , i> = -<t? + , (46) 

and allows to write oj in terms of it 

u = -(<P + fa) . (47) 
Finally, the third equation of (H) gives the celebrated Hirota equation [24] 

. + <h - 0i - , fi. v . <p l + 0x + (j>\ + <p 

sm = tan — tan — sm . 48) 

2 2 2 2 y J 
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